Development of a deformation-tunable quadrupolar microcavity 
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We have developed a technique for realizing a two-dimensional quadrupolar microcavity with 
its deformation variable from 0% to 20% continuously. We employed a microjet ejected from a 
noncircular orifice in order to generate a stationary column with modulated quadrupolar deformation 
in its cross section. Wavelength red shifts of low-order cavity modes due to shape deformation 
were measured and were found to be in good agreement with the wave calculation for the same 
deformation, indicating the observed deformation is quadrupolar in nature. 

PACS numbers: 05.45.Mt, 42.55.Sa 



I. INTRODUCTION 



Unlike circular microcavities showing isotropic output 
emission, asymmetric resonant cavities (ARC's)P, |2 Q 
provide directional output along with fairly high cavity 
quality factors up to certain degree of deformation. In 
addition, the internal ray dynamics of ARC's exhibits an 
evolution from regularity to chaos as the degree of defor- 
mation increases. With these properties, ARC's serve not 
only as a useful resonators for optoelectronics devices but 
also as a versatile tool for studying ray and wave chaos 
i. 

In arc's, the internal ray dynamics strongly depend 
on the degree of deformation, and so do mode proper- 
ties such as output directionality |^ and cavity quality 
factor. Whispering gallery modes and their ultrahigh 
cavity qualtiy Q factors of circular microcavities can be 
associated with regular ray dynamics in the cavities. As 
the deformation is increased gradually, this dynamical 
regularity begins to be broken following Kolmogorov- 
Arnold-Moser(KAM) scenario in phase space. The emis- 
sion pattern also experiences a dramatic change, from 
isotropic emission into coUimated tangential emission at 
high-curvature points on the surface. When the deforma- 
tion is further increased, even KAM tori are broken and 
the internal ray dynamics becomes more chaotic. How- 
ever, this ray picture is not fully applicable to typical 
microcavities with their radial size ranging from several 
microns to tens of microns. In this region, it is not well 
understood how the chaotic ray dynamics is related to 
both actual modes of cavity and output emission pat- 
terns of these modes. 

If one can tune the deformation of a microcavity from 
a perfect circle to a final deformed form continuously. 
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one can follow the evolution of modes and thus can eas- 
ily identify the origin and mode numbers of observed 
modes of ARC. Such identification allows us to perform 
direct comparison with wave calculation and help us to 
obtain better understanding on the connection between 
the mode distribution, cavity quality factors and out- 
put directionality of these modes and the aforementioned 
chaotic ray dynamics. 

Notwithstanding the advantage of electrical puirmin^ 
capability, microcavities made of semiconductor [T 

are rigid under room temperature and thus cannot 
be changed in shape continuously. In order to have var- 
ied shapes, separate cavities with different deformation 
should be employed. However, it is almost impossible to 
make connections among observed modes from the sepa- 
rate cavities. In addition, the cavity quality factors are 
limited to 10^ ~ 10^ without cumbersome annealing pro- 
cess even for fairly large microcavities with diameters of 
hundreds of microns. For such large cavity sizes the free 
spectral range of modes are comparable to or larger than 
the linewidth of modes and thus individual modes are 
difficult to be resolved. 

On the other hand, surface-tension-induced microcav- 
ities (STIM's), such as liquid dro plet s llo| . Hquid 
micro i et fllJ . siHca micro spheres [Tj, [Tjl and micro- 
toroid [ij can provide ultrahigh-Q factors at relatively 
small sizes due to their surface smoothness. STIM's 
have been employed in studying shape-deformation ef- 
fects 0illll|l3- Among various STIM's, liquid microjet 
technique is particularly suitable for generating deformed 
two-dimensional microcavities. 

In the present work, we have realized a stable station- 
ary microjet ARC with continuously variable deforma- 
tion by a dapting the well-developed droplet generation 
technique [lall7lll8| and precisely controlling the hydro- 
dynamics of the microjet. Our liquid microcavity offers 
another benefit that the concentration of gain molecules 
for lasing and other spectroscopic measurements can be 
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easily adjusted. Although it is only subject to optical 
pumping, the liquid microjet ARC is nevertheless a pow- 
erful system for studying the aforementioned ray-wave 
correspondence or the quantum chaos in the microcav- 
ity. 

This paper is organized in the following way. In Sec. 
II, we describe a method to realize deformed microcavi- 
ties with a set of discrete degrees of deformation by us- 
ing a stationary surface-modulated microjet column. A 
method to obtain quadrupole cavities with continuously 
tunable deformation by controlling the pressure of the 
jet is discussed in Sec. III. Spectral analysis of cavity 
modes in deformed microcavities is presented in Sec. IV, 
where we demonstrate that the wavelength red shift due 
to shape deformation only can be extracted from raw ex- 
perimental data. The observed shifts vs. the measured 
degrees of deformation are shown to be in good agree- 
ment with the results of wave calculation. 



II. METHOD TO OBTAIN DISCRETE DEGREE 
OF DEFORMATION 

An ARC can be obtained from a horizontal cross sec- 
tion of a microjet, which is made of ethanol (71=1. 361) 
doped with Rhodamine B dye as gain molecules. A sim- 
ilar microjet generator with a circular orifice was intro- 
duced previously 0. However, in the present study the 
microjet employs a noncircular orifice as a nozzle. 

We have reported elsewhere that such a noncircu- 
lar orifice induces a microjet to form a stationary tidal 
column as depicted in Fig. ^a). Viscosity of the liquid 
makes the oscillation amplitude of the column rapidly 
damped along the jet propagation direction. The shape 
of the jet column can be approximated in the cylindrical 
coordinates by the following time-independent equation 



9,z) 



l + r]Q exp 



27r 



z + M cos 26* 



(1) 

where a is the mean radius and rja is a seed deforma- 
tion parameter, t is the decay rate, T is the period of 
oscillation and ^ is an initial phase of oscillation. Both 
the decay rate and the period of oscillation are approxi- 
mately constant depending on the properties of the liquid 
and the orifice size. The jet velocity Vz is assumed to be 
uniform across the entire jet [20l |. 

A two-dimensional ARC, specifically a quadrupole- 
deformed microcavity (QDM), is obtained by selecting 
cross-sectional planes located at the extrema z„'s (n=-l, 
0, 1, 2,. . . in Fig.^a)) of the amplitude oscillation, where 
the cavity boundary is given by 



r(0,z„) = a[l + (-)"?7„cos20] 



where 
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FIG. 1: (a) Model for the deformed microjet column, (b) 
Fabrication procedure of a noncircular orifice. The side view 
shows that the inner walls are inclined. The bottom view is 
a real image of a noncircular orifice. 



It is our convention to denote the cross-sectional plane 
or anti-nodal plane at z„ as Dn (n = 1,2,3...). Ac- 
cording to our convention, Dl corresponds to the second 
minimum of amplitude oscillation from the nozzle in Fig. 
^a). Equation ^ describes a quadrupole of deforma- 
tion 77 = 77,1, and thus we can obtain QDM's at Dl, D2, 
D3, etc. The major axis of QDM's at Dl, D3, D5, etc is 
rotated by 90° with respect to that of QDM's at D2, D4, 
D6, etc. 

Degree of possible three-dimensional (3D) effect due to 
a finite thickness of the region to be used in experiment 
around the cross-sectional plane, which is about 10 mi- 
crons, can be estimated in the following way. Under our 
experimental conditions to be described below, a typical 
value of VzT is 280 microns, which is much larger than 
the mean radius a of 15 microns. The possible variation 
in the cavity size over the 10 /zm region in the z direction 
is then estimated to be less than 0.1% or 0.02 micron 
in diameter, and therefore, the 3D effect can be safely 
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neglected. 

The degree of deformation of QDM's can be deter- 
mined from their diffraction patterns made by an inci- 
dent probe laser as described in Refs. ^.nd "T^. De- 
tails of diffraction pattern measurement and data anal- 
ysis are found in Ref. The degrees of deformation 
of QDM's are found to be 32% (Dl), 19.2% (D2), 13.3% 
(D3), 10.2% (D4), 6.5% (D5), etc under 1.8 bar of ejec- 
tion pressure. 

Let us now describe the fabricating process of the non- 
circular orifice, a key component in our system. Figure 
n^b) depicts the fabrication procedure. A tip of pyrex 
tube is placed vertically in the center of a Nichrome heat- 
ing coil. As the tube melts down, the bottom part of the 
tube droops. At that moment, the inner walls of the 
tube are inclined, resembling the letter 'V as depicted 
in the side view of Fig. ^b). While being heated, the 
both sides of the tip are pressed with a tweezer. After 
slightly grinding off and polishing the tip we can obtain a 
noncircular orifice on it. The diameter can be controlled 
from several microns to hundreds of microns by amount 
of polishing. 



III. METHOD TO OBTAIN CONTINUOUSLY 
TUNED DEGREE OF DEFORMATION 

According to Eq. (1), the microjet boundary at the 
nozzle (z — 0) is given by 



r{0, 0) = a (1 + 770 sin ^ cos 29) 



(4) 



which should match the boundary of the orifice. In 
fact, the boundary of the orifice usually contains higher- 
order components of small magnitude. These higher- 
order components have larger damping rates than that 
of the quadrupole, and therefore, by the time of the jet 
reaches Dl location, the microjet surface becomes well 
approximated by Eq. 

Since the degree of deformation at the nozzle is more 
or less prescribed by the shape of the orifice, rji = rjo sin ^ 
can be regarded to be fixed. However, r]o, called a seed 
deformation, can still be changed under this constraint, 
and therefore, the deformation 77,1 at anti-nodal plane 
Dn {n = 1,2,3...), proportional to 770 as in Eq. ((JJ, 
can be modified accordingly. The seed deformation rjo 
can be tuned by the ejection pressure of the jet. For 
example, with an ejection pressure increased to 2.0 bar, 
we obtained a new set of QDM's with deformations of 
35% (Dl), 22% (D2), 14.4% (D3), 11% (D4), 7.6% (D5), 
etc, increased from 32% (Dl), 19.2% (D2), 13.3% (D3), 
10.2% (D4), 6.5% (D5), etc at 1.8 bar. 

The pinching of the nozzle in one direction in the fab- 
rication process has introduced different inner wall slope 
and this differential wall slope in turn induces different 
initial radial velocities, which makes it possible to tune 
the cavity deformation. The initial radial velocity just 
outside the orifice can be obtained from Eq. J^l by a 




FIG. 2: Cross sectional view of the deformation oscillation 
along the z axis, given by Eq. (0 with ^ = 0. Initial deforma- 
tion 770 sin^) at z = is predetermined by the orifice. As 
the jet velocity is increased, the initial phase ^ changes in 
such a way that cos ^ is proportional to the velocity (see Eq. 
Q). Consequently, the deformations at anti-nodal planes, 
Dl, D2, D3 . . ., can be increased as the jet velocity. The pa- 
rameter values used for these plots are T — 28/^s, r = 39^s, 
a = 15/im, K=0.1 and 77i = 0.7. 



taking time derivative: 



VriO) 



/27ra 



rjo cos £, cos 29 



(5) 



where we have neglected the damping effect in taking the 
time derivative. Due to the guided compression by the in- 
ner wall of the nozzle, the radial velocity inside the nozzle 
is proportional to the ejection velocity. Therefore, we can 
assume that the initial radial velocity of the quadrupole 
oscillation is also proportional to the jet velocity in the 
z direction just outside the orifice: 



.(61) = Kv^ cos 26* 



(6) 



where is a constant depending only on the inner geom- 
etry of the nozzle and the properties of the liquid. It is 
independent of the ejection pressure. From Eqs. Q and 
© , we can see that cos ^ must be proportional to Vz '■ 



rjQ cos^ = KvzT /2'Ka 



(7) 



With this and the relation 77^ = 770 sin ^, we then obtain 
the following result showing that the seed deformation 
can be controlled by the jet ejection velocity. 




(8) 



For the representative values of the parameters in Eq. 
0, r ~ 30 /iS, Vz ~ 10 m/s, a ~ 15 /im, rji ~ 0.7 and 
K ~ 0.1, the ratio {2TTar]i/ KvzT) ~ 2.2. The jet ejection 
velocity is easily controlled by the ejection pressure. 
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It is noted that the oscillation of the quadrupole defor- 
mation along the z direction is analogous to a damped 
harmonic oscillator with a non-zero launching velocity. 
The dependence of ^ and consequently the amplitude of 
deformation oscillation on the jet velocity Vz is illustrated 
in Fig. 2. In addition, the contribution due to the damp- 
ing in the initial radial velocity is much smaller than that 
of the harmonic oscillation. For the parameter values in 
Fig. 2, the ratio is T/27rr - 0.1. 

Deformation tuning by the jet ejection pressure is 
experimentally demonstrated in Fig. |21 where the de- 
grees of deformation r\n at several anti-nodal planes Dn's 
(n = 2,3,4,5) are plotted as a function of the ejection 
pressure P, which is varied from 1.4 bar to 2.6 bar by a 
0.2-bar step. A stable microjet is obtained in this pres- 
sure range. By choosing a proper combination of the 
anti-nodal planes Dn and the jet pressure P, we can tune 
the deformation to any value from 5% (D5, 1.4 bar) to 
28% (D2, 2.6 bar) as shown in Fig. 3. Some deformation 
can be achieved with more than one combination of Y)n 
and P. 

Small deformation from 0% to 4% can be obtained 
around high anti-nodal planes (06^^010) under a fixed 
pressure. The small deformation can be more precisely 
measured with spectroscopic methods in the next sec- 
tion than with the diffraction patterns since the diffrac- 
tion patterns for such slightly deformed QDM's are not 
distinct enough to be used for determining their defor- 
mations reliably. 



IV. SHAPE CONFIRMATION BY SPECTRUM 
ANALYSIS 

A. Method and Principle 

Low-order whispering gallery modes in a circular mi- 
crocavity reside very close to the cavity boundary. As the 
microcavity is deformed from a circle to an asymmetric 
shape while the internal area is preserved, these low-order 
modes experience wavelength shifts. Since the low-order 
modes still reside very close to the cavity boundary even 
in a deformed cavity, these wavelength shifts can be re- 
lated to the change in the perimeter of the cavity. 

The wavelength shifts due to the change in the perime- 
ter of a great circle of a micro droplet under volume- 
preserving shape deformation have been investigated be- 
fore [2lll23| . In that case, the observed wavelength shift 
was mostly due to the size increase of the great circle 
under consideration and thus did not give definitive in- 
formation on the shape deformation itself for the great 
circle. In the present work, however, we consider shape 
deformation for a two-dimensional microcavity with its 
internal area preserved, and therefore, by measuring the 
wavelength shift, we can find out how much degree of 
deformation the cavity has. 

Spectra of cavity modes can be observed in several 
different ways. We use the cavity modified fluorescence 
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FIG. 3: (a) Variation of deformation parameter according 
to the ejection pressure of the jet at D2, D3, D4 and D5 
positions, (b) For deformation below 4%, the deformation is 
tunable by selecting different anti-nodal planes at a fixed jet 
pressure of 1.4 bar. 



of cavity medium. Specifically, we measure the fluores- 
cence of dye molecules resolved in the microjet. Due to 
the cavity quantum electrodynamics effect, the fluores- 
cence spectrum exhibits enhancement wherever the cav- 
ity modes are located and thus it replicates the resonance 
lines of the cavity modes. 

The setup for spectrum measurement is similar to that 
of Ref. [ill. The microjet is excited by an argon-ion 
pump laser at 514 m. The pump laser is focused with a 
cylindrical lens into a thin profile with a thickness of 10 
micron in the z direction so that a thin slab can be se- 
lected to be a two-dimensional microcavity in the micro- 
jet. The fluorescence emitted from this region is collected 
by an objective lens and delivered to a spectrometer with 
a charge-coupled-device detector. The observed fluores- 
cence spectrum shows several sequences of peaks with 



5 



(a) 



3 

c 

U 
0} 



^ , 111. 



where 



2.6 

2.4 
2.2 
2.0 
1.8 
1.6 
1.4 



590 600 610 620 630 640 
Wavelength (nm) 

(b) 



15 
12 
9 
6 
3 

-3 
-6 
-9 
■12 



4-++- 



1.4 1.6 1.8 2.0 2.2 2.4 2.6 
Pressure (bar) 



FIG. 4: (a) Cavity-modified fiuorescence spectra observed at 
D3 as the jet pressure is varied. A particular mode of a good 
visibility marked by arrows is followed as the jet pressure is 
changed, (b) Observed wavelength shifts of the particular 
mode in (a) , measured for the jet pressure varied at a small 
interval, are denoted by filled squares. Wavelength blue shifts 
due to the area contraction as the pressure is increased are 
represented by filled circles. The wavelength shifts due to 
the deformation only, denoted by crosses, are obtained by 
subtracting the latter from the former. Error bars are smaller 
than the symbol sizes. 



well-defined free spectral ranges (FSR's) . The modes sep- 
arated by the same FSR are the modes of the same mode 
order. In this way, we can identify the mode orders of all 
of the observed modes. 

A low-order mode, which has a good visibility, can 
be tracked in the spectrum as the deformation is var- 
ied. Some examples are shown in Fig. E^a). Since the 
low-order mode reside very close to the cavity bound- 
ary as discussed above, the resonance wavelengths of 
these modes are approximately proportional to the cavity 
perimeter S, given by 



S = 



r{9) = 



V1 + W2 



(1 + 7? cos 26*) , 



(10) 



(9) 



and r] is the quadrupole deformation parameter, i.e., one 
of 77„'s at anti- nodal planes. The square root factor in the 
denominator is a normalization factor needed to make the 
area of the quadrupole conserved regardless of rj. A solid 
curve in Fig. shows the red shift for given deformation 
according to Eqs. (9) and (10). This phenomenon will be 
refered to as 'pure deformation effect' from now on. 



B. Experimental results 

Figure 0Ja) shows the observed spectra for various jet 
pressure values. The peaks marked by arrows correspond 
to the same mode undergoing wavelength shifts. Anti- 
nodal plane D3 was chosen and the peak was tracked 
from 1.4 bar to 2.6 bar. The deformation parameters 
measured by the diffraction technique were 11.1%, 12.0%, 
13.3%, 14.4%, 16.9%, 17.7% and 19.4%, respectively, in 
the ascending order of the listed pressure values. In Fig. 
4(b), the wavelength shifts with respect to the wave- 
length at 1.4 bar, measured at 0.25 bar interval, are 
plotted in filled squares. Contrary to our expectation, 
they show blue shifts as the deformation increases from 
1.4 bar to 2.0 bar, beyond which they start to show red 
shifts. We have found that the blue shift occurs because 
the cross-sectional area contracts as the ejection pressure, 
and accordingly the ejection velocity, increases due to the 
Bernoulli effect. 

We checked how much area contraction occurs by mea- 
suring spectra at anti-nodal plane D8. Since the defor- 
mation there is less than 2%, the deformation effect can 
be neglected, and thus the wavelength shift of a partic- 
ular mode due to area contraction can be measured as 
a function of the ejection pressure Since area contrac- 
tion accompanies reduction in perimeter, it leads to a 
blue shift. Filled circles in Fig. EJb) represent the ob- 
served wavelength shifts at D8 as a function of the ejec- 
tion pressure. They indeed exhibit pure blue shifts. The 
area contraction is believed to be completed immediately 
after the jet leaves the orifice. 

By subtracting the wavelength shifts due to the area 
contraction from the overall wavelength shifts, we obtain 
the wavelength shifts due to the deformation effect only. 
The results are shown in Fig. 01b) as crosses, exhibiting 
pure red shifts as expected. 

The range of deformation obtained from D3 is limited, 
from 11% to 19.4%. Smaller deformation can be obtained 
from D4 and D5. Since the measurement of wavelength 
shift requires that the same mode be followed as the de- 
formation changes, we need to find the same mode in D4 
and D5 as the one used in D3. We can directly compare 
the spectrum obtained from D3 at 1.4 bar with the spec- 
trum from D4 at 2.0 bar since both give 11% deformation 
as shown in Fig. Ola). By using the fact that the same 
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mode should have the same FSR, we can identify the 
mode in the spectrum from D4 that is the same as the 
mode tracked in the spectrum from D3. Once the same 
mode is found, we can adjust the pressure and extend the 
wavelength shift measurement from 2.6 bar down to 1.4 
bar for D4. Similarly, we can connect the measurement 
from D4 to D5. In this way, we can follow one particu- 
lar mode with a good visibility across several anti-nodal 
planes and measure the wavelength shift for deformation 
from 5%(at 1.4 bar in D5) to 19.4% (at 2.6 bar in D3). 



C. Small deformation case 

Deformation smaller than 5% can be obtained from 
the region of jet from D5 and above. Since the deforma- 
tion is so small there, the surface profile is nearly flat, 
and thus any thin slice in this region can serve as a two- 
dimensional microcavity. We can easily track one mode 
in the observed spectra while we scan the excitation lo- 
cation continuously with the ejection pressure of the jet 
fixed. The observed wavelength shifts for various excita- 
tion location on the jet at 1.4 bar are shown in Fig.|5fa). 
The oscillatory behavior is due to the oscillation of the 
cavity perimeter according to the oscillation of the degree 
of deformation as in Eq. (^). The overall underlying red 
shift is due to the slight area expansion as the jet goes 
up. The vertical velocity of the jet Vz decreases slightly 
as the jet goes up due to air friction and the gravity. 
Since the flux of the jet should be conserved, the cross 
sectional area slightly expands as the velocity drops. 

The red shift can be easily extracted from the observed 
shifts. The local minima in the observed shifts must cor- 
respond to no deformation (i.e., circle) and they are fit 
well by a straight line as shown in Fig.l^b). These local 
minimum values are the red shifts due to the area expan- 
sion. The amount of red shift over a distance from one 
anti-nodal plane to the next is only 0.3 nm. After sub- 
tracting the amount of red shifts corresponding to this 
linear fit from the observed shifts, we obtain the shifts 
due to the deformation only as shown in Fig. Efc) . 

The observed oscillatory behavior conforms to an ef- 
fective deformation parameter rj^s given by 



??cff(^:) = ??oexp 



27r 

sm I -—z + S^ 

VzT 



(11) 



The zero shift occur when the sine factor vanishes (no 
deformation) whereas the local maxima occurs when the 
sine factor is ±1 (at the anti-nodal planes). The exper- 
imental points for the degree of deformation below 5% 
in Fig. 1^1 are obtained from the peak shift values at the 
anti- nodal planes DB'-^^DIO. 

We can calculate the perimeter of a quadrupole defined 
by Eq. (10) and from this perimeter we can calculate 
wavelength shift due to deformation only with respect to 
the wavelength for 77 = 0. The result is shown as a solid 
curve in Fig. El In our experiment, we can measure the 
degree of deformation as a function of the jet pressure 




1.2 1.4 
Position (mm) 

FIG. 5: (a) Wavelength shifts in the small deformation re- 
gion. Oscillating pattern shows the deformation effect and 
the underlying slope is due to the area expansion as the jet 
goes up. (b) Wavelength shifts due to the area expansion can 
be obtained by joining the local minima in (a), corresponding 
to no deformation, (c) Wavelength shifts due to deformation 
only is obtained by subtracting the shifts due to the area ex- 
pansion in (b) from the total shifts in (a). Solid lines are the 
fits. The fitting parameters in (c) are ?7o=0.7, r = 35 /.ts, 
T — 28 fis, ^=2.1 for Vz=10 m/s. Error bars are smaller than 
the point sizes. 



by using the diffraction technique as mentioned before. 
We can also measure the wavelength red shift due to 
deformation only for a given jet pressure as discussed 
above. We then obtain a series of pairs of a red shift and a 
degree of deformation for various jet pressure values. For 
small deformation, we measure the red shift for a fixed 
pressure but at different position on the jet as discussed 
above. The results are shown as filled squares in Fig. 
The experimental results are in good agreement with the 
calculation assuming quadrupole deformation. 

Small deviation observed for large deformation is due 
to the fact that considering the low order modes followed 
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in spectra as a ray circulating very close to the cavity 
boundary is just too much simplification. We can solve 
the wave equation for a given quadrupole cavity with the 
same size parameter as in the experiment and find res- 
onance wavelengths of the same mode observed in the 
experiment. Details of mode calculations are given in 
Ref. '23^. The resulting wavelength shifts from the wave 
calculations are shown in Fig. 6 as open triangles, which 
also show about the same amount of deviation from the 
solid curve by the simple model. The observed wave- 
length shifts match well with the results of the wave cal- 
culation, and therefore, we conclude that the shape of our 
microcavity can be regarded as a quadrupolar deformed 
microcavity within the limit of our experimental error. 



FIG. 6: Observed wavelength red shifts vs. the degree of de- 
formation measured by the diffraction technique. Solid line is 
the calculated red shift from the perimeter of a quadrupole, 
open triangles are the result of wave calculation for the same 
cavity and filled squares are the experimental results. Dot- 
ted line is the calculated red shift for an ellipse with its area 
conserved and with the same deformation as the quadrupole. 
Vertical error bars are smaller than the point size. 
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